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ABSTRACT 

We show that the observed scatter in the relations between the mass of supermassive black holes 
(SMBHs), Mbh, and the velocity dispersion Csph or mass Mgph of their host spheroid, place interesting 
constraints on the process that regulates SMBH growth in galaxies. When combined with the observed 
properties of early-type SDSS galaxies, the observed intrinsic scatters imply that SMBH growth is reg- 
ulated by the spheroid velocity dispersion rather than its mass. The Mbh-Mgph relation is therefore a 
by-product of a more fundamental Mbh-Csph relation. We construct a theoretical model for the scatter 
among baryon modified dark matter halo profiles, out of which we generate a population of spheroid 
hosts and show that these naturally lead to a relation between effective radius and velocity dispersion of 
the form i?sph ~ 6kpc((Jsph/200kms~^)^-^ with a scatter of '^0.2dex, in agreement with the correspond- 
ing projection of the fundamental plane for early type galaxies in SDSS. At the redshift of formation, 
our model predicts the minimum scatter that SMBHs can have at fixed velocity dispersion or spheroid 
mass under different formation scenarios. We also estimate the additional scatter that is introduced 
into these relations through collisionless mergers of purely stellar spheroids at z < 1 . We find that the 
observed scatter in the Afbh-Csph and Mbh-Afgph relations preclude the properties of dark matter halos 
from being the governing factor in SMBH growth. The apparent relation between halo and SMBH mass 
is merely a reflection of the fact that massive halos tend to host massive stellar spheroids (albeit with a 
large scatter owing to the variance in formation histories) . Finally, we show that SMBH growth governed 
by the properties of the host spheroid can lead to the observed values of scatter in the Mbh-cTsph and 
Mbh-Afsph relations, only if the SMBH growth is limited by momentum or energy feedback over the 
dynamical time of the host spheroid. 

Subject headings: cosmology: theory - galaxies: formation 



1. INTRODUCTION 

Supermassive black holes (SMBHs) are a ubiquitous 
constituent of spheroids in nearby galaxies (e.g. Kor- 
mandy & Richstone 1995). A decade ago it became ap- 
parent that the masses of SMBHs correlate with the lumi- 
nosity of the host spheroid (Kormandy & Richstone 1995). 
Subsequently, other correlations with substantially less in- 
trinsic scatter have been discovered. These include cor- 
relations between the SMBH mass (Afbh), and the mass, 
Msph (Magorrian et al. 1998; Haering & Rix 2004), steUar 
velocity dispersion, ctsph (Ferrarese & Merritt 2000; Geb- 
hardt et al. 2000), and concentration (Graham et al. 2002) 
of its host spheroid. The tightest relation, with intrin- 
sic scatter 6 ^ 0.275 ± 0.05dex, appears to be between 
SMBH mass and velocity dispersion (Tremaine et al. 2002; 
Wyithe 2005). 

These relations must hold important clues about the 
astrophysics that regulates the growth of a SMBH and 
its impact on galaxy formation. While much attention 
was dedicated towards interpreting the power-law slope 
of the Mbh-Csph relation via a slew of analytic, semi- 
analytic and numerical attempts to reproduce it (e.g. Silk 
& Rees 1998; Wyithe & Loeb 2003; King 2003; Miralda- 
Escude & KoUmeier 2005; Adams et al. 2003; Sazonov et 
al. 2005; Di Matteo et al. 2005), little attention has been 
directed towards interpreting the constraints that its in- 
trinsic scatter might place on our understanding of SMBH 
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growth (Robertson et al. 2005). Moreover, agreement be- 
tween data and theory is a necessary but not sufficient 
condition. A model that reproduces the observations is 
not necessarily unique. The various successful attempts 
to model the quasar luminosity function assuming differ- 
ent physical models attest to this fact. Our goals in this 
paper are to use the observed scatter in the Mbh~crsph and 
Mbh-Afsph relations in addition to their power-law slope as 
a diagnostic of SMBH formation physics, and to constrain 
a range of possible models of SMBH formation. 

Throughout the paper we adopt the set of cosmologi- 
cal parameters determined by the Wilkinson Microwave 
Anisotropy Probe (WMAP, Spergel et al. 2003), namely 
mass density parameters of rim = 0.27 in matter, flf, = 
0.044 in baryons, n\ = 0.73 in a cosmological constant, 
and a Hubble constant of Hq = 71 kms~^ Mpc""'^. 

2. INTRINSIC SCATTER IN THE OBSERVED mbh-CTsph AND 
"^bh""^sph RELATIONS 

Tremaine et al. (2002) have compiled a list of spheroids 
with reliable determinations of both SMBH mass and cen- 
tral velocity dispersion (defined as the luminosity weighted 
dispersion in a slit aperture of half length i?sph)- These 
SMBHs show a tight correlation between Afbh and agph 
(Gebhardt et al. 2000; Ferrarese & Merritt 2000). Re- 
cently, Wyithe (2005) found that the Mbh-Cgph relation 
within the sample of Tremaine et al. (2002) shows ev- 
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idence for a deviation from a pure power-law behavior, 
and obtained an a-posteriori probability distribution for 
the scatter in the relation of 6 = 0.28 ± 0.05. The best-fit 
log-quadratic Mth-Csph takes the form 

logio(Mbh) = a -h /Jlogio (cTsph/200km/s) 



+ p2 [logio (CTsph/200km/s)]' 



(1) 



where /3 — 4.2 ±0.37 corresponds to the traditional power- 
law relation and P2 = 1.6 ± 1.3 quantifies the departure of 
the local SMBH sample from a pure power-law. 

A correlation has also been found between the mass of 
the spheroid component, Mgph, and the SMBH (Magorrian 
et al. 1998). The sample compiled by Haering & Rix (2004) 
(which overlaps predominantly with the Tremaine et 
al. 2002 sample) was studied by Wyithe (2005) who found 
that like the Afbh fsph relation, the Mbh^Msph relation 
may depart from a uniform power-law, and that the scat- 
ter in the Mbh-Afgph relation is larger than in the Mbh- 
(Tsph relation by 50% ((5sph = 0.41 ± 0.07). The best-fit 
log-quadratic Mbh-Mgph takes the form 

logio(Mbh) = CKsph + /Ssphlogio (Msph/IO^^M©) 

+ /32,sph [logio (Msph/lO^Me)] ' , (2) 

where /3sph = 1.15 ± 0.18 and /Ja.sph = 1-12 ± 0.14. 

In this paper we use the values of intrinsic scatter {5 
and (5sph) in the local SMBH sample to place constraints 
on the process that regulates SMBH growth in galaxies. 

3. WHICH RELATION IS MORE FUNDAMENTAL: mbh-Csph 
OR TObh-TOsph? 

The SMBH mass is observed to correlate tightly with 
both (Tsph (Tremaine et al. 2002) and Mgph (Haering & 
Rix 2004), implying that SMBH growth is regulated by 
properties of the spheroid. It is therefore natural to ask 
which parameter among these two is responsible for setting 
the SMBH mass? We address this question empirically 
without specifying a mechanism for the SMBH growth. 

Let lis suppose first that the A'/bh A'-^sph relation is the 
fundamental one. There would still be a relation between 
Mbh and crgph since Msph and agph are related. Bernardi et 
al. (2003) have compiled values of -Rgph, and o-gph for the 
sample of early-type galaxies in the SDSS. We can gen- 
erate a sub-sample of galaxies within a narrow range of 
(Tsph and find the dynamical mass Mgph = ^ph-^sph/G = 
2fTgpjji?sph/G for each, where Vgph is the circular velocity at 
iisph and we have assumed an isotropic velocity dispersion 
c^sph = Vsph/V^- We can then use the observed relation 
Mbh oc Msph with its intrinsic scatter of Sgph = 0.41 ±0.07, 
to find the corresponding scatter in the Mbh^fsph relation. 
The resulting distributions of residuals in the SMBH mass 
are plotted on the left panel of Figure 1, together with the 
1-sigma range in the observed scatter (grey region) of the 
-^bh-CTsph relation. If the Mbh-Mgph relation were funda- 
mental, then the scatter in the Mbh^Cgph relation would 
have been 5 = 0.46dex, well in excess of the observed value 
5 = 0.275 ±0.05. This implies that the Mbh-M^ph relation 
is not fundamental. 

On the other hand if we suppose that it is the Mbh- 
fjsph relation that is fundamental, then there would still 
be a relation between Mbh and Mgph. We generate a sub- 
sample of galaxies from Bernardi et al. (2003) with Mgph = 



2c7gpjji?sph/G in a narrow range near Mgph = lO^^M©. We 
can then use the observed Mbh-fsph relation (with intrin- 
sic scatter 5 = 0.275±0.05) to find the corresponding scat- 
ter in the A'/bh Mgph relation. The resulting distributions 
of residuals in SMBH mass are plotted in the right hand 
panel of Figure 1, together with the 1-sigma range in the 
observed scatter (grey region) of the Mbh-Mgph relation. 
If the A^bh-Cgph relation is fundamental, then the scat- 
ter in the Mbh-Mgph relation should be ^sph — 0.40dex, in 
good agreement with the observed value ^gph = 0.41±0.07. 
This implies that the Mbh-Cgph relation is more fundamen- 
tal than the Mbh-Mgph relation, with the latter being an 
incidental consequence of the correlation between asph and 
M,ph. 

4. INTRINSIC SCATTER IN THE mbh-CTsph RELATION AND 
MODELS FOR SMBH EVOLUTION 

In the previous section we showed that the observed 
scatter in the Mbh fsph and Afbh Msph relations, when 
combined with the scatter among spheroid properties im- 
plies that the Mbh^Csph relation is fundamental for SMBH 
growth, with the Mbh-Mgph relation being incidental due 
to correlations between spheroid parameters. This in turn 
implies that if we have a model for the properties of the 
host spheroid, then we can deduce the mode of SMBH 
growth by comparing the calculated scatter in the modeled 
A^bh o'sph and Afbh Afgph relations with observations. In 
this section we constrain the astrophysics of SMBH growth 
by computing minimum values for the scatter in the Mbh- 
CTsph relation using various models for the regulation of 
SMBH growth. 

4.1. Rotation curves, Adiabatic cooling and the 
Fundamental Plane 

Our goal is to model the scatter in the A'/bh fsph and 
Mbh-Mgph relations. To accomplish this goal we must have 
a representative model of the host spheroids as well as the 
scatter in their parameters. This is achieved by comput- 
ing the rotation curve that results from cooling of baryons 
inside a dark-matter halo. 

The virial radius iiyir and velocity Vvir for a halo of mass 
Afhaio at redshift z are 
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where = [l + {nA/nm){l + z)~^]-\ = 187r2±82d- 
39d^ and d = — 1 (see Barkana & Loeb 2001 for more 
details). The relation of the circular velocity at the virial 
radius to the velocity dispersion at smaller galactic radii 
requires specification of the mass density profile. In this 
work we adopt the Navarro, Frenk & White (1997, here- 
after NFW) profile for the dark matter. In addition to V^ir 
and Mvir the NFW profile is defined by the concentration 
parameter c, which is the ratio between the virial radius 
and a characteristic break radius, c = r^ir/rg- The me- 
dian value of c depends on halo mass Mhaio and redshift 
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iogio(Mbh/Mo) iogio(^^bh/^^o) 

Fic;. 1. Left: The predicted intrinsic cumulative distribution of residuals (logj^Q(Mijh/Mo), where Mq is the mass corresponding to 
the mean local relation) in M^h at constant Csph = 200 km s^^. In this case the observed M^h— Mspjj relation with intrinsic scatter 
of <5sph = 0-41 ± 0.07dex was assumed to be fundamental. Right: The predicted intrinsic residual scatter in the local Mt,h at constant 
Mjph = 10^^ Mq assuming the observed Mth-Csph relation with intrinsic scatter of 5 = 0.275 ± O.OSdex to be fundamental. In each panel the 
range of the observed scatter (1-sigma) is shown as the shaded grey region for comparison. Zero scatter would have been represented by a 
step function. 




Fig. 2. — Scatter plot of Rgph vs Usph from cooled profiles inside dark matter halos (black-dots). For comparison, the SDSS sample (Bernardi 
et al. 2003) is shown by the grey dots. The thick solid and dashed lines show the mean observed relation Rgph = 6(crsph/200)^'^ and the level 
of observed scatter (±0.2dex). The thin solid line shows a linear relation. 



z. Based on N-body simulations, Bullock et al. (2001) (see 
also Wechsler et al. 2002) have found a median relation 

for c, 

with a scatter of Alogc = 0.14 dex. 

The parameters describing the spheroid are its charac- 
teristic radius i?sph and velocity dispersion o-gph- In this 
paper wc assume that the gas available to cool in the halo 
makes most of the mass within the effective radius of the 
spheroid Mgph with a density distribution described by a 
Hernquist (1990) profile. The cooling of baryons modifies 
the density profile of a galaxy and hence its rotation curve. 
In particular, the velocity dispersion in the luminous por- 
tion of the galaxy is substantially larger than would be 
expected from an NFW profile. Gnedin et al. (2004) have 
studied the contraction of an NFW halo due to baryon 
cooling in a cosmological context. They find that tradi- 
tional adiabatic contraction does not provide a good fit. 
However they introduce an alternative modified contrac- 
tion model and provide fitting formulae for contracted pro- 



files as a function of halo mass, concentration parameter, 

final characteristic radius for the baryons, and the cooled 
baryon fraction. From this contracted profile we can com- 
pute the velocity dispersion fTgph at the characteristic ra- 
dius i?sph- However the cooled profile, and hence the value 
of CTsph obtained from the cooled profile depends on the 
value adopted for i?sph- Defining to be the fraction of 
the total galaxy mass that makes the spheroid (including 
the cooled baryons), we can break this degeneracy by iden- 
tifying the spheroid mass mdMhaio with the effective virial 
mass Mgph through 2a^^^^Rsph/G = M^ph = rridMhaio- 
With this second relation we are able to solve uniquely for 
i?sph and CTsph within a specified dark matter halo (with 
parameters Mhaio c, z, and ma). We choose a probability 
distribution that is flat in the logarithm of md over the 
range 0.015 < md < 0.15. The upper value corresponds to 
the case where all baryons in the halo cool to form the mass 
of the spheroid [so that Mgph = (f^b/f^m)Afhaio]j and the 
range under consideration spans values smaller by up to 
an order of magnitude relative to this extreme case. Using 
this formalism we generate a sample of model spheroids. 
This sample is compared first to the observed fundamental 



4 



plane of early- type galaxies (below), and then used to dis- 
cuss the scatter in models for the Mbh-Csph and Mbh-Afgph 
relations (§ 4.2). 

Figure 2 shows a scatter plot of i?sph vs dsph for cooled 
profiles with randomly selected parameters. The distribu- 
tion of formation redshift is dictated by the merger trees 
of halos and corresponding fate of pre-existing galaxies 
within these halos (Eisenstein & Loeb 1996). It is known 
empirically that star formation is rather minimal between 
< z < I for elliptical galaxies (Bcrnardi ct al. 2003). 
Since we find that the relation between i?sph and crgph is 
not very sensitive to the precise probability distribution 
of collapse redshifts, we adopt in this figure a distribu- 
tion of collapse redshifts that is flat between 1 < z < 3. 
The resulting points show a correlation with signific;ant 
scatter. For comparison, the relation between the velocity 
dispersion cTsph of an early type stellar system and its ef- 
fective radius for galaxies from the SDSS (Bernardi et 
al. 2003) is shown by the grey dots. A fit to this observed 
correlation has the form 

which is plotted as the thick solid line in Figure 2. The 
scatter about the observed relation, ±0.2dex, is bracketed 
by the pair of dashed lines. Also plotted to guide the eye 
is a linear relation (thin line). On comparison with the 
results of Bernardi et al. (2003), we see that the formalism 
reproduces the observed behavior Rgph oc (jj^- Moreover 
the size of the predicted scatter about the mean relation 
is ~ O.lTdex, in good agreement with the observed value 
of ~ 0.2dex. This result is not very sensitive to the (un- 
known) distribution chosen for rrid- The model does not 
contain any free parameters. Nevertheless, in addition to 
the scatter an power-law slope of the i?sph^crsph relation, 
our prescription also predicts its normalization. 

For completeness we note that the relation described by 
equation (6) depends only on dynamical properties and 
not on the details of star-formation. However elliptical 
galaxies follow a 3-parameter fundamental plane, with the 
scatter around the median i^sph-CTgph relation parameter- 
ized by surface brightness rather than a 2 parameter rela- 
tion (Djorgovski & Davis 1987). The scatter aroimd the 
fundamental plane is substantially smaller than around 
the i?sph~crsph relation shown in Figure 2. Bernardi et 
al (2003) find that the fundamental plane has the form 

n ^ ^1.49±0.05 r-0.75±0.01 /7^ 

where the surface brightness is Igph oc L/R-^^^^, and show 
that this requires the mass-to-light T to have a depen- 

1 /3 

dence on iJgph of the form F oc R^!^^- A complete model of 
the fundamental plane would need to explain this depen- 
dence of the mass-to-light ratio, in addition to the relation 
^sph oc CTgp^. However we do not expect SMBH growth to 
depend on F, and do not attempt to model this (orthogo- 
nal) parameter. 

Finally, we note that the non-linearity of the relation be- 
tween i?sph oc CTgpfj is surprising since it is derived within 
the context of CDM dark-matter halos for which the rela- 
tion between the virial radius and velocity is i?vir oc K'ir at 
any given redshift. In the following section we discuss the 
growth of SMBHs inside host spheroids in light of the ob- 
served scatters in the Mbh-Csph and Mbh-Mgph relations. 



4.2. Models for SMBH Evolution 

Three observations of the relation between SMBHs and 
their hosts have motivated di£Fc;rent classes of models to 
describe the growth and evolution of SMBHs through ac- 
cretion. First, the observation of the Magorrian et al. 
(1998) relation that the mass of the SMBH follows the 
mass of the spheroid has motivated models where the mass 
of the SMBH accretes a constant fraction of the available 
gas following a major merger (e.g. Haiman & Loeb 1998; 
Kauffmann & Haehnelt 2000). We refer to this scenario 
as case-I. Second, the observation of the Afbh cTsph relation 
implies that the SMBH growth is regulated by the depth of 
the gravitational potential well of its host spheroid. We re- 
fer to this scenario as case-II. Third, there is evidence from 
the observation of the Mbh^o-gph relation in quasars, that 
the Mbir o"sph relation does not vary with redshift (Shields 
et al. 2003). This motivates a revision of case-II, to include 
regulation of the SMBH growth over the dynamical time of 
the system (e.g. Silk & Rees 1998; Haehnelt, Natarajan & 
Rees 1998; Wyithe & Loeb 2003). We refer to this as case- 
III. On the other hand if niomcntmn rather than energy is 
conserved in the transfer of energy in a quasars outflow to 
the cold galactic gas, then rather than a SMBH regulated 
by binding energy, we have a SMBH mass regulated by 
the binding energy divided by the virial velocity (Silk & 
Rees 1998; King 2003; Begelman 2004; Murray, Quataert, 
& Thompson 2005). In analogy with cases-II and III, the 
SMBH mass may be regulated by the total momentum of 
the surrounding gas, or by the total momentum divided by 
the systems dynamical time. We refer to these as cases-IV 
and V, respectively. 

We therefore test five hypotheses for each of the Mbh- 
fTsph and Mbh-Mgph relations. In the previous section we 
described a model that reproduces the observed behavior 
of i?sph oc CTgp^, with a scatter of ~ 0.2dex. The agreement 
with the observed projection of the fundamental plane 
(Bernardi et al. 2003), gives us confidence that the model 
provides a sufficiently accurate framework within which 
wc can discuss the scatter in the Mbh-crgph and Mbh~-^sph 
relations. Below we list the details of all cases under con- 
sideration: 

• Case-I: the mass of the SMBH saturates at a con- 
stant fraction of the mass of the spheroid, Mbh oc 
Msp\, oc tr^pj^i^sph. 

• Case-II: the mass of the black-hole grows in pro- 
portion to the binding energy of baryons in the 
spheroid. Taking a constant fraction of the cold-gas 
mass to be material that must be expelled from the 
spheroid during the feedback we find that the black- 
hole mass is therefore Mbh oc Msp^al^^ oc a^^^Rsph- 

• Case-Ill: the mass of the SMBH is determined by 
the mass for which accretion at the Eddington limit 

provides a constant fraction of the binding energy 
of the baryons in the spheroid over a constant frac- 
tion of the spheroid's dynamical time. Thus, the 
black- hole mass scales as Mbh oc -E'b/(-Rsph/o'sph) oc 
A^sph^-^phZ-Rsph oc (jfph. 
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Fig. 3. — Predicted residual scatter (logjQ(Mbh/Mo), where Mo is the SMBH mass corresponding to the mean relation) relative to the 
best-fit log-quadratic relations. Results are shown as a function of cr^pij (upper row) and Mgpij (lower row) for models where the SMBH 
growth is regulated by properties of the spheroid. The dashed grey lines show the level of observed scatter. 
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Fig. 4. — Predicted intrinsic scatter for models where the SMBH growth is regulated by properties of the spheroid. In each case we show 
the predicted cumulative distribution of residuals (logj^oCA^bh/Aft)! where Mo is the SMBH mass corresponding to the mean relation) at 
constant asph = 200km/s (left) and constant A/sph = IO^^Mq (right). The dashed, dot-dashed, solid, dotted, and dot-dot-dot-daslied lines 
correspond to cases I, H, HI, IV and V respectively. In the left panel the lines for cases-Ill and V, and for cases-I, II and IV are coincident. 
In each panel the range of the observed scatter (1-sigma) is shown as the shaded grey region for comparison. 



• Casc-IV: As in case-II, but with the momentum 
rather than the energy of the outflow couphng to 
the gas in the spheroid, yielding Mbh oc MsphO-gph oc 

^sph-^sph- 

• Case-V: As in case-Ill, but with the momentum 
rather than the energy of the outflow coupling to 
the gas in the spheroid over a dynamical time, yield- 
ing MsphCT^pj^/iisph oc CT^^ph- 

Studies of the local SMBH inventory suggest that most 
of mass in SMBHs was accreted during a luminous quasar 
phase (e.g. Yu & Tremaine 2002; Shankar et al. 2004), 
with a potentially significant contribution from an addi- 
tional dust-obscured accretion phase (Martinez-Sansigre 
et al. 2005). If the fraction of obscured quasars is inde- 
pendent of redshift, then the quasar luminosity function 
(Fan et al. 2004; Boyle et al. 2000) can be used as a proxy 



for the distribution of SMBH formation redshifts. Using 
this redshift distribution and the formalism outlined in the 
previous sub-section, we can estimate the slope and scat- 
ter in the Mbh-fsph and M^h-Msph relations under the five 
different scenarios for the regulation of SMBH growth. 

For each of the five cases we perform Monte-Carlo sim- 
ulations of SMBH growth. The resulting residuals in the 
.^bh-Csph and Mhh^Msph relations are plotted in Figure 3 
relative to the best fit log-quadratic relations (equations 1- 
2) as a function of crgph (upper row) and Mgph (lower row) 
respectively. For each of the five cases, wc also compute 
distributions of residuals (in dex) relative to the mean 
SMBH mass at constant asph and constant Mgph- The 
resulting distributions are plotted in Figure 4. The values 
of scatter are listed in the 1st and 2nd columns of Table 1. 
The ratio of the scatter at constant M^ph and at constant 
cTsph are listed in column 3. The power-law slopes /? and 
Psph in each case are listed in columns 4 and 5. 
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Table 1 

Predicted scatter and logarithmic slope in the Mbh-o-gph relation and the Mbh-Msph relation. 

spheroid halo 





5 






P 


/^sph 


5 


^sph 


Ssph/S 


/5 


Psph 


Case-I 


0.17 (0.275) 


0.0 (0.22) 


(0.8) 


3.4 


1.0 


0.24 (0.275) 


0.29 (0.31) 


1.2 (1.1) 


3.3 


1.0 


Case-II 


0.17 (0.275) 


0.10 (0.24) 


0.6 (0.87) 


5.3 


1.6 


0.36 


0.48 


1.3 


5.4 


1.6 


Case-Ill 


0.0 (0.275) 


0.25 (0.37) 


oo (1.35) 


5.0 


1.5 


0.37 


0.53 


1.4 


5.5 


1.5 


Case-IV 


0.17 (0.275) 


0.05 (0.22) 


0.29 (0.8) 


4.2 


1.3 


0.39 


0.38 


1.0 


4.3 


1.3 


Case-V 


0.0 (0.275) 


0.20 (0.34) 


oo (1.24) 


4.0 


1.2 


0.29 


0.42 


1.4 


4.4 


1.3 



Prom the above discussion we see that under the as- 
sumption of a unique coupling efficiency between the 
quasar output and the surrounding spheroid, cases-Ill and 
V imply a perfect Mbh^cTsph relation. Cases-I, II and IV 
do not produce perfect Mbh o"sph relations even under this 
unique coupling assumption due to the scatter in the i?sph- 
iTsph relation. Similarly, the scatter in the i?sph ^sph rela- 
tion leads to scatter in the Afbh^-^sph relation for cases-Ill 
and V. As a result, cases-Ill and V predict a scatter in the 
Mbh-Afsph relation that is larger than in the Mbh-Csph 
relation, in agreement with observations. 

Obviously in these five cases we have neglected many 
additional sources of scatter. These include, but arc not 
limited to, the fraction of the Eddington limit at which 
the SMBH shines during its luminous phase, the efficiency 
of coupling of feedback energy or momentum to the gas in 
the host spheroid, and the fraction of the system's charac- 
teristic size for which the dynamical time should be com- 
puted. In all five cases the minimum values of the scatter 
are smaller than the observed 5 = 0.275. However in cases 
I, II and IV, the models predict a scatter in the Mbh^Afgph 
relation that is smaller than in the Mbh-CTgph relation while 
the observations show the opposite. On the other hand if 
SMBHs are regulated by the binding energy or momentum 
of gas in the spheroid per dynamical time of the spheroid 
(cases III & V), then the minimum scatter in the Mgph- 
(Tsph relation is reduced to zero. Cases-Ill and V therefore 
predic;t that the scatter in the Mbh-Mgph relation should 
be larger than in the Mbh-Csph relation as observed. 

The scatter in the projection of the fundamental plane 
onto the i?sph o"sph plane therefore allows us to differenti- 
ate between SMBH growth that is regulated by the mass 
(Case-I), binding energy (Case-II) or momentum of gas 
(Case-IV) in the spheroid, and SMBH growth that is reg- 
ulated by energy or momentum feedback over a dynamical 
time of the spheroid (Cases-Ill and V). In the latter cases 
the scatter in the Mbh-Mgph relation is increased by the 
large scatter in the spheroid radius, Rgph- On the other 
hand in Cases-Ill and V, i?sph cancels out in the division 
of mass by dynamical time in the determination of Afbh at 
constant agph- 

If the predicted value for the minimum scatter in the 
-^bh-f sph relation is smaller than the observed value of 
S = 0.275 dex, then there is room in the model for ad- 
ditional random scatter to account for varying Eddington 
ratio, outflow geometry, dust composition, and other fac- 
tors. In each case we therefore add random scatter in the 
formation process at a level which results in the predicted 
scatter in the Mbh-Csph relation being equal to the ob- 



served value oi 6 = 0.275 dex. This value, the correspond- 
ing value predicted for the scatter in the Mbh-Mgph rela- 
tion ((5sph), and the corresponding ratio (Ssph/S) are listed 
in parentheses in Table 1 . While cases-I, II and IV predict 
ratios of scatter between the Mbh^Msph and Mbh-Csph re- 
lation that are smaller than unity, case-III predicts a ratio 
of Ssph/S ~ 1-4, and case-V a ratio of 6sph/S ~ 1.2, in good 
agreement with the observed value (Ssph/S = 1.5 ± 0.35). 
We therefore conclude that SMBH growth is likely regu- 
lated by feedback over the spheroid's dynamical time. 

A further possible discriminant between models is the 
power-law slope of the Mbh-CTsph relations. Cases HI & 
V, which satisfy constraints from the observed scatter in 
the Mbh-Csph and Mbh-Afsjph relations, produce power-law 
-^bh-f sph relations, with SMBH mass in proportion to the 
velocity dispersion raised to the fifth and fourth power re- 
spectively. This could be compared with the power-law 
value from Tremaine et al. (2002) oi (3 — 4±0.3 for galaxies 
with CTsph ~ 200 km s~^, a comparison which at first sight 
this appears to support case-V. However Wyithe (2005) 
has found evidence for a power-law slope that varies with 
(Jsph from /3 ~ 4 near (Tsph ~ 200 km s^^, to /? ^ 5 near 
fsph ~ 350 km s~^ (see Eq. 1). The slope of the Mbh^ 
Mgph relation is observed to be close to unity (Haering & 
Rix 2004). Of the cases (HI and V) that produce an ac- 
ceptably small scatter for the Mbh Csph relation, we find 
that case-III yields /?sph ~ 1-5, while case-V leads to a 
value of Psph ~ 1-2. However Wyithe (2005) also finds ev- 
idence for a varying power-law in the Mbh-Afgph relation, 
though not at high significance. 

Since the models described in this paper each predict 
power-law relations between Mbh and tisph, the residuals 
with respect to the log-quadratic fit therefore show cur- 
vature as a function of Csph- Cases-Ill and V both show 
a power-law slope that agrees with the best-fit relation at 
some values of (Tgph, which renders discrimination between 
models based on their predicted power-law slope difficult. 
It is beyond the scope of this paper to discuss detailed 
models of the Mbh^Cgph relation; however the observed 
curvature may point to a turn-over from momentum to 
energy conservation in the reaction between the outflow 
and surrounding gas. Alternatively, there could be a ve- 
locity dependent efficiency of feedback, which could change 
the slope of the relation. In summary, based on the ob- 
served values of (5 = 0.275 ± 0.05, Ssph = 0.41 ± 0.07, only 
models corrcisponding to Case-III or Case-V are accept- 
able. Examples of this class include the models of Silk & 
Rees (1998); Wyithe & Loeb (2003); King (2003), di Mat- 
teo et al. (2005), or Murray et al. (2005). 
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4.3. Additional scatter from dissipationless mergers after 
the quasar epoch 

Our model computes the minimum scatter in the Afbh 
iJsph relation during a time when the SMBH can grow via 
gas accretion. At late times a SMBH may find itself in 
an environment where there is no remaining cold gas. In 
this regime a merger of two galaxies will proceed via col- 
lisionless dynamics. The two SMBHs may coalesce once 
they enter the gravity-wave dominated regime, but since 
there is no cold gas the SMBH will not grow during a 
quasar phase. It would require detailed numerical simu- 
lations to discover whether or not the Mbh-Cgph relation 
holds following a collisionless merger. However, unlike the 
situation where feedback is still at work, in a collisionless 
merger the stellar spheroid is not sensitive to the existence 
of the SMBH (at least at radii comparable to -Rgph)- More- 
over the total SMBH mass is nearly conserved as long as 
SMBH binaries coalesce due to interactions with the sur- 
rounding stars and the emission of gravitational radiation 
(Begelman, Blandford, & Rees 1980). The growth of the 
SMBH and the properties of the spheroid following a col- 
lisionless merger should be independent in the sense of the 
-^bh-fsph relation. Thus, we would expect collisionless 
mergers to increase the scatter in the Mbh^Csph relation, 
leaving our estimates of the minimum scatter intact. 

N-body simulations of the behavior of stars as collision- 
less particles have been performed by Gao et al. (2004) who 
found that the inner mass density profile (interior to some 
fixed physical radius) is unaffected by mergers, implying 
that the velocity dispersion interior to any given radius 
remains the same after a merger. When two spheroids 
merge, their combined stars cover a larger radius (because 
their mass increases while the inner mass profile remains 
unchanged), and this changes the value of crgph. We may 
therefore use the fact that the inner profile remains in- 
variant in order to estimate the scaling between i?sph and 
CTsph in the regime where purely collisionless mergers oc- 
cur. Based on the simulations of Gao et al. (2004), we 
assume a universal NFW mass profile for the total mass 
(dark matter-l-stars) irrespective of the merger history and 
find what mergers would do to the Mbh-fgph relation when 
the total mass in stars and in SMBHs is conserved. If the 
inner density profile maintains the NFW shape of p cx l/r, 

then cTgpjj oc r oc M^^^, i.e. Mgph oc af^^, similar to the 
Faber- Jackson (1976) projection of the fundamental plane 
of spheroids. 

Collisionless mergers will change the average normaliza- 
tion of the Mbh-CTsph relation. To see why, suppose that we 
have Np equal mass progenitors at redshift z, each with ve- 
locity dispersion CTgph.p and spheroid mass Mgph.p. At red- 
shift z, the SMBHs obey Mbh.p = Cpcfph pi where /3p is the 
slope and Cp is a constant. The final SMBH and spheroid 
masses at = are Mbh,f = A^pAfbh,p and Mgph.f = 
A^pAfsph.p respectively. Using the above scaling the final 

1 /4 

velocity dispersion is (Jgph.f = o"sph,p^p • We therefore 
find Mbh.f = iVpMbh.p = iVpCpaJ-, p = CpTVp^-'^^/VJ', 
Thus the normalization of the Mbh^cTsph relation at fixed 

CTgph is changed by a factor Np ^^^^ through collision- 
less mergers. Note that the slope of the Mbh-crgph relation 
/3p is preserved if A^p is independent of o-gph- However 



the number of progenitors is a function of halo mass and 
redshift, and so the change in normalization could be a 
function of o-gph- 

In addition to changing the normalization of the Mbh- 
(Jsph relation, collisionless mergers will also introduce scat- 
ter. In the limit of a perfect correlation between Mgph and 
CTsph, and where the number of progenitors (A^p) is con- 
stant, the argument in the previous paragraph shows that 
collisionless mergers would lead to an Mbh-o-sph relation 
with no additional scatter beyond that introduced at the 
formation redshift. However the scatter in the i^sph-fsph 
correlation combined with the relation Mgph oc (J^p^Rsph 
implies that there is ^ 0.2dex of scatter in Mgph at fixed 
(Jsph- Moreover, different galaxies have different merger 
histories and therefore a different number of progenitors. 
Scatter among the properties of the initial building blocks 
at redshift z therefore leads to scatter in the local Mbh- 
Cgph relation even if the Afbh cTsph relation at z were per- 
fect. This scatter introduced through collisionless mergers 
will therefore add scatter to the local Mbh-fgph relation 
beyond that intrinsic to the formation process itself. 

To ascertain the quantitative effect of collisionless merg- 
ers on scatter in the Mbh o"sph relation, we generate merger 
trees of dark-matter halos using the method described in 
Vollonteri, Haardt & Madau (2003). Based on the merger 
trees we find the A'^p progenitor halos at z ^ 1.5 that lead 
to a halo of known mass at z ~ 0. Using the formalism 
outlined in § 4.1 we then determine the values of crgph, 
i?sph and Mgph for the spheroids populating these progen- 
itor halos. We also populate the spheroids with SMBHs 
of mass Afbh according to the perfect Afbh~crsph relations 
that arise from cases-Ill and V. The final SMBH mass re- 
siding in the halo at z = is Afbh,f = X^i=o -^bh,i- It 
is embedded in a spheroid of mass Afsph.f — X]i=o -^sph.i- 
Based on the above scaling for the inner p cx l/r den- 
sity profile of an NFW halo we may estimate the value 
of velocity dispersion corresponding to the final spheroid, 

Csph,!' = o-sph,o(Afsph,f/Afsph,o)^^'', where Msph,o and (Tsph^o 
are the mass and velocity dispersion of the largest pro- 
genitor. In Figure 5 we show the scatter introduced into 
a perfect Afbh-o'sph relation originating at z = 0.5 (left), 
z = 1.5 (center) and z = 2.5 (right) by collisionless merg- 
ers between those redshifts and z = 0. As discussed above 
this scatter arises as a result of scatter in the relation be- 
tween CTsph and Afgph in the progenitors. The upper panels 
show results for Case-Ill (/3 = 5), while the lower pan- 
els show results for Case-V (/3 = 4). The scatter intro- 
duced is roughly independent of Usph and takes values of 
5 ~ O.ldex, 5 ^ 0.2dex and 5 ^ O.Sdex for mergers origi- 
nating atz = 0.5, z = 1.5 and z = 2.5, respectively. Thus 
galaxies that become devoid of gas at higher redshift lead 
to a larger scatter in the Afbh Csph relation, because these 
galaxies undergo more collisionless mergers by z = than 
a galaxy which becomes devoid of gas only at late times. 

The stars that populate massive galaxies appear to be 
older than those in low mass galaxies (Kauffmann et al. 
2003). The cold gas reservoir that made these stars must 
have been depleted at a higher redshift for the progenitors 
of high-CTsph galaxies. We might therefore expect more 
scatter in the Afbh~crsph relation at large CTgph- For equal 
mass mergers we have shown that the normalization of 
the Afbh-cTgph relation at a fixed agph changes by a factor 
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Fig. 5. — Scatter plots of the predicted residuals (logjg(Mbh/Mo), where Mq is the SMBH mass corresponding to the mean power-law 
relation at rcdshift z) in that result from coUisionless galaxy mergers. Each point represents a different realization of a coUisionless 

merger tree beginning with a perfect power-law Mbij-o"sph correlation at z = 0.5 (left panel), z = 1.5 (center panel) and z = 2.5 (right panel). 
We assume Case-Ill /3 = 5 (upper row) and Case-V (3 = 4 (lower row). 



~ Np Thus, for /3p = 5 (Case-Ill), we find that the 

amplitude of the Mbh-CTsph relation should be reduced by 
coUisionless mergers, while for /3p = 4 (Case-V) the am- 
plitude should be preserved. This behavior is seen in Fig- 
ure 5. Moreover, more massive galaxies undergo a larger 
rate of major mergers. Figure 5 shows that in Case-V the 
change in normalization of the Mbh-Csph relation through 
coUisionless mergers is more significant in high mass than 
in low mass galaxies, as expected. CoUisionless mergers 
could therefore lead to a reduction in the steepness of the 
observed Mbh-Csph relation if /3p > 4. 

In summary, in order to satisfy the constraint that the 
-^bh-Csph relation have a scatter at ^; = that is smaller 
than ~ 0.3dex, the intrinsic scatter in the relation at the 
formation rcdshift should be smaller than ^0.2 dcx, so as 
to allow for the additional scatter introduced through col- 
lisionless mergers. Cases-Ill and V meet this requirement. 

4.4. Redshift dependence of the Mbh-fsph o,nd Mbh-Msph 

relations 

Recent evidence suggests that the Mbh-tTsph relation is 
preserved out to high redshift (Shields et al. 2003), while 
SMBHs make a larger fraction of their host spheroid mass 
at higher rcdshift (Rix ct al. 2001; Groom ct al. 2004; 
Walter et al. 2004). This behavior is reproduced in mod- 
els with scenarios of the form case-Ill or case-V. Figure 6 
shows a scatter plot of the predicted residuals in Mbh 
vs z at constant asph = 200 km s~^ (upper row) and 
Msph — 1O"'^^M0 (lower row) for each of the five models. 
The residuals are normalized relative to the mean relation 
at 2: = 3 in each case. While there is no evolution in the 
Mbh-CTsph relation for cases-Ill and V, we see that SMBHs 
are predicted to be an order of magnitude more massive 
with respect to their host spheroid at z ~ 6 than they 
are at 2; 1 in agreement with observations. In contrast, 
models for case-I, II and IV predict that the SMBH mass 



should decrease by an order of magnitude between z = 
and z = 6 at constant tJsph, while not evolving significantly 
at constant -Mgph. The observed evolution of SMBH mass 
with redshift therefore supports cases III and V as the sce- 
nario for SMBH growth. This result supports the findings 
of § 4.2 based on the scatter in local relations. 

4.5. Do dark-matter hales play a role in SMBH 

evolution? 

Attempts to reproduce the observed luminosity func- 
tion of quasars associate the mass of the SMBH with the 
properties of the host dark-matter halo. This paradigm 
allows the abundance of SMBHs to be traced either 
in semi-analytic or mmicrical models (e.g. Haiman & 
Locb 1998; Hachnolt, Natarajan & Rocs 1998; Wyithe & 
Loeb 2003; VoUonteri, Haardt & Madau 2002). Indeed, 
Ferrarese (2001) has inferred a relation between the masses 
of SMBHs and their host dark-matter halo. Is it possi- 
ble that the halo rather than the spheroid regulates SMBH 
growth ? 

Since we have computed spheroid properties within a 
specified dark matter halo, we are in a position to discuss 
the role of the dark matter halo in regulating the SMBH 
growth. In addition to the five cases listed above for SMBH 
growth within a spheroid, we also try five analogous cases 
for the formation of SMBHs governed by dark matter halo 
properties. For each of the additional five cases, we com- 
pute the distribution of residuals (in dex) relative to the 
mean relation as a function of (Tsph and -Mgph via a Monte- 
Carlo algorithm for SMBH formation, and calculate the 
variance at each of constant CFsph = 200 km s~^ and con- 
stant Msph = 10"'^"'^ Mq. Below we list the details of each 
case. The resulting distributions of residuals are plotted in 
figure 7. The values of scatter are listed in the 6th and 7th 
columns of Table 1. The ratio of the scatter at constant 
Msph and at constant agph are listed in column 8. Values 
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Fig. 6. — Scatter plots of the predicted residuals (log]^Q(Mi,h/A^o), where Mq is the SMBH mass corresponding to the mean relation) in 
Mbh as a function of z. Results axe shown at constant (Tgph = 200 km s~^ (upper row) and constant M^ph = 10^^ (lower row) for models 
where SMBH formation is governed by spheroid properties. In each of the five cases shown the residuals are relative to the mean relation at 
z = 3. 
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Fig. 7. — Predicted scatter of residuals (logjQ (Afbh/A<fo) where Mo is the SMBH mass corresponding to the mean relation) relative to 
the best-fit log-quadratic relations. Results are shown as a function of cTsph (upper row) and Mgph (lower row) for models where the SMBH 
growth is regulated by properties of the dark matter halo. The dashed grey lines show the level of observed scatter. 



for the slopes /? and /?sph are listed in columns 9 and 10. 

• Casc-I: the mass of the SMBH forms from a con- 
stant fraction of the baryonic component of the halo 
mass, Mbh oc (Ob/nm)Mhaio- 

• Case-II: the mass of the black-hole grows in propor- 
tion to the binding energy of baryons in the halo. 
For an NFW profile with a concentration parameter 

c we get Mbh cx (rJb/Om)MhaioKk.fc, where 

^ c l-l/(l + c)^-21n(l + c)/(l+c) 

2 [c/(l + c) - ln(l + c)]2 ■ ^ ' 

• Case-Ill: the mass of the SMBH is determined by 
the mass for which accretion at the Eddington limit 
provides the binding energy of baryons in the halo 
over a constant fraction of the halo's dynamical 



time (Wyithe & Loeb 2003). We therefore have 

M, , ( /Hm ) Mhalo Vyi, /c f T/5 

iWbh ^X H^^^/V^ir ^Jcl'vir- 

• Case-IV: same as case-II, but with the momen- 
tum rather than the energy of the outflow cou- 
pling to the gas in the spheroid. We then find 

Mbh oc (f2b/^^m)MhaloKir/c 

• Case-V: same as case-IV, but with the momentum 
rather than the energy of the outfiow coupling to 
the gas in the spheroid over a dynamical time. We 
then find Mbh oc oc /cKi- 

In these five cases we have again neglected many possi- 
ble causes of scatter. However in case-I, where the SMBH 
growth is regulated by halo properties, the minimum value 
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Fig. 8. — The predicted distribution of intrinsic residuals (logiQ(Mi,h/A^o) where Mq is the SMBH mass corresponding to the mean relation) 
in the local M^h at M^aio = 10^^ Mq. The M^h-Afhaio relation has larger intrinsic scatter than the Mbh-Csph relation. 



of the scatter in the Mbh-Csph relation is slightly smaller 
than the observed 5 = 0.275. We have added random 
scatter to the model for case-I in order to bring the pre- 
dicted scatter in the Mbh-CTgph up to the observed value. 
This value and the corresponding prediction for Sgph and 
S/Ssph are listed in parentheses. Case-I predicts a scatter 
in the Mbh Msph relation that is similar to the Mbh-cTsph 
relation while the observations suggest the latter to be sig- 
nificantly smaller. While case-1 cannot be ruled out based 
only on the predicted minimum scatter in the Mbh-Csph 
relation at the formation redshift, the small allowance for 
additional expected scatter from the aforementioned astro- 
physical sources renders it unlikely, particularly when the 
additional scatter of 0.1 — O.Sdex from collisionless mergers 
at low redshift is accounted for. In cases II, III, IV and 
V where the SMBH growth is regulated by halo proper- 
ties, the minimum values of the scatter in the Mbh^crgph 
relation are larger than the observed (5 = 0.275. 

While scatter within models of spheroid regulated 
SMBH growth is not sensitive to the distribution of rud, 
the predicted scatter in models of halo regulated SMBH 
growth decreases as the assumed range of decreases. 
However, in order to reduce the predicted scatter below the 
^ 0.2dex threshold at the formation redshift (which would 
allow for additional scatter to be introduced through col- 
lisionless mergers), we find that the allowed range around 
Wd = 0.05 would need to be smaller than ±0.025, which is 
implausibly narrow. We therefore conclude that it is the 
spheroid rather than the dark-matter halo which drives 
the evolution of SMBH mass. 

4.6. The Mbh^Afhaio relation 

We have demonstrated that the tight relation between 
Mbh and fJsph implies that it is the spheroid rather than 
the halo which governs the growth of SMBHs. However it 
is clear that since there is an Mbh-o-sph relation, and since 
larger halos will, on average, host bulges with larger cen- 
tral velocity dispersions, there should also be a correlation 
between SMBH and halo mass. Ferrarese (2001) has found 
such a relation. Since it is not possible to measure the dark 
matter halo mass directly, halo masses for galaxies in the 
local sample were inferred via a maximum circular veloc- 
ity estimated from o-gph based on an empirical relation. 
It is therefore difficult to estimate the scatter in the Mbh- 



A^haio relation observationally. Here we predict the scatter 
in the Mbh-Mhaio relation at a fixed value of Mhaio- We 
compute the distribution of residuals via a Monte-Carlo 
method as before. We choose Mhaio = IO^^Mq and find 
the distribution of values for fJsph, and hence the distri- 
bution of Mbh using the observed Mbh-Csph relation. The 
resulting distribution is plotted in Figure 8. The variance 
is Shaio — 0.4 dex. Thus the tightness of the Mbh o'sph rela- 
tion suggests that the Mbh-Mhaio correlation is incidental 
to the fundamental relation between the SMBH and its 
host spheroid. Note that this variance is computed at the 
time of the SMBH formation. The surrounding dark mat- 
ter halo could continue to grow after the supply of cold gas 
to the SMBH had ceased. This is consistent with our con- 
clusion that SMBHs grow in proportion to the properties 
of the spheroid rather than the halo. Indeed one finds mas- 
sive dark-matter halos in X-ray clusters, which must have 
increased their velocity dispersion well beyond the corre- 
sponding SMBH growth (due to the lack of cooling flows 
in cluster cores). This late-time growth of dark-matter 
halos increases considerably the scatter in the Mbh-Mhaio 
relation. 

5. CONCLUSION 

We have investigated the implications of intrinsic scat- 
ter in the local relations involving SMBHs for models of 

SMBH formation. Using the sample of spheroid proper- 
ties from SDSS (Bernardi et al. 2003) we first examined 
empirically the fundamental parameter describing SMBH 
growth. The observed scatter in the Mbh-o-gph relation 
is (5 = 0.275 ± 0.05, while the Mbh-Mgph relation has a 
larger observed scatter of Ssph = 0.41 ± 0.07. Assuming 
that the Mbh^Mgph relation is fundamental, we use the 
SDSS spheroid sample to compute the resulting scatter in 
the Mbh-cTsph relation. We find that this procedure results 
in a scatter of the Mbh-^sph relation that is too large to 
be reconciled with observation. Alternatively, one might 
assume that SMBH growth is determined by CTsph rather 
than Mgph. In this case we used the SMBH sample to 
compute the resulting scatter in the Mbh-Mgph relation, 
and found agreement with the observed scatter. We there- 
fore conclude that SMBH growth is governed by agph, and 
that the observed correlation between Mbh and Mgph is a 
by-product of the relation between Mgph and cTsph- 
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Theoretical models for SMBH forraation must repro- 
duce several observational constraints: (i) the scatter in 
the local Mbh-crgph relation is 5 = 0.275 ± 0.05 dex, im- 
plying that at the time of formation the scatter should be 
smaller than ^ 0.2dex to allow for additional scatter in- 
troduced by coUisionless mergers of galaxies since z ~ 1 or 
earlier; (ii) the scatter in the Mbh-Mgph relation is larger 
than in the Mbh-f sph relation (this result is maintained as 
additional scatter from coUisionless mergers is introduced 
after SMBH formation); and (in) The Mbh o"sph relation 
is preserved out to high redshift. We find that these con- 
straints are only met by models where SMBH growth is 
regulated by feedback on the gas feeding the SMBH over 
the spheroid dynamical time. Other models lead to scatter 
in the Mbh" o"sph relation that are too large or scatter in 
the Mbh-Afgph relation that is smaller than the Mbh-Csph 
relation. In addition, other models lead to a SMBH mass 
that drops with increasing redshift at a fixed velocity dis- 
persion. The feedback in successful models can be either 



in the form of energy or momentum transfer between the 
quasar and the galactic gas, leading to power-law slopes 
in the Mbh-Csph relation of /3 = 4 or /3 = 5, respectively. 
Both of these slopes are permitted by the local sample 
(Wyithc 2005). 

The above constraints do not permit SMBH growth to 
be governed by the properties of the dark- matter halos. 
Such models lead to scatter in the Mbh-o'sph relation that 
are too large and/or scatter in the Mbh Mgph relation that 
is smaller than the Mbh-o'sph relation. The relation be- 
tween Mbh and the halo mass (Ferrarese 2001) has a large 
scatter (~ 0.4dex) and is most likely a by-product of the 
correlation between halo mass and Csph. 
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